Free energy analysis of system comprising biased atomic force 
microscope tip, water meniscus and dielectric surface 
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We are concerned with free energy analysis of the system comprising an AFM tip, water meniscus, 
and polymer film. Under applied electrostatic potential, the minimum in free energy is at a dis- 
tance greater than the initial tip-substrate separation in the absence of potential. This equilibrium 
distance, to, mostly depends on the tip bias V and cantilever spring constant k s , where as varia- 
tions of io is less pronounced with respect to the dielectric constants, and polymer film thickness. 
Polarization of water meniscus under the AFM tip appears to be the dominant factor enabling the 
creation of mechanical work for tip retraction. 

PACS numbers: 68.37.Ps (Atomic force microscopy), 68.47.Mn (Polymer surfaces) 
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I. INTRODUCTION 

The goal of this work is to analyze the free energy of 
the system comprising an electrically biased AFM tip, 
water meniscus, and dielectric surface. It is important 
to identify the key factors leading to the effective tip- 
substrate repulsion through an equilibrium free-energy 
balance. Complete understanding of the process is com- 
plicated by the fact that the specific spatial details of 
the tip-surface contact profile as well as cantilever mo- 
tion with applied bias is exceedingly difficult to observe 
and quantify; and thus is not well understood. Although 
in this work we concentrate on the theoretical consid- 
eration of this process for polymer films in the context 
of the AFM-assisted electrostatic nanolithographyji the 
equilibrium free energy balance and associated conclu- 
sions are generally applicable to other material systems 
including semiconductors, bio-macromolecules and self- 
assembled monolayers. 



II. MODEL DESCRIPTION 

Consider the schematic representation in Figure ^ as 
summarizing the key components involved in the process. 
The system comprises the AFM tip, a dielectric film, and 
the water meniscus filling the junction. The tip is pre- 
sented as a sphere of radius R attached to the cantilever 
with the spring constant k s . The region of interest is 
near the tip apex. Tip-surface separation t confines this 
region in the vertical direction, and the cylinder of the 
radius R c defines its lateral dimensions. The total free 
energy of the system can be presented as the sum of four 
terms: 




F = F 1 +F 2 + U + F tl 



(1) 



Here, F\ and F 2 are Helmholtz's free energy contribu- 
tions related to polarization of the dielectric film, and 



FIG. 1: Conceptual presentation of AFM tip-polymer surface 
junction. 



water meniscus respectively; U is the energy of the con- 
tracted spring, and F t8 is the energy of electrostatic at- 
traction between the tip and the film. The free energy 
of the system is chosen to be zero for an unbiased tip at 
t = and thus F represents the change in free energy. 
Free energy densities variations at the tip/water and the 
water/polymer interfaces due to electric field, and possi- 
ble charge carriers' generation are not accounted in the 
Eq. (TJ>. 

We consider the terms F\ and F 2 separately, since the 
variation of dielectric medium entropy with respect to ex- 
ternal electric field is greater for water than that for poly- 
mer. Additionally, we neglect the term related to Van 
der Waals interactions between the tip and the surface: 
Ui/dw — —AR/6t because the typical value of Hamaker 
constant for condensed matter is in the order of 10" 19 J 2 
resulting in energy less than 10~ 29 J, which is found to 
be negligible in comparison to the other free energy terms 



2 



in JTJ. 

The dielectric film is considered as a non-conductive 
polymer with a macroscopic dielectric constant t\. The 
first term thus can be written as: Fy = f fidV, where 

V 

the free energy density fy is 

h = (2) 

where i? = —V<f> is the external electric field, eo is the 
dielectric permittivity in free space, and we assume that 
ei is independent of E even for the substantial fields (~ 
10 9 Vm -1 ) considered. 

To account for the polarization of a liquid medium, the 
free energy of water has a form 



Fo = -k B T\nZ 



(3) 



where Z is the total partition function, ks is Boltz- 
man's constant, and T is the temperature. The parti- 
tion function can be written for the canonical ensemble 
of N molecules of liquid as Z = Z c Q N /N\A Here, Q 
is the partition function of a molecule, and Z c is the 
configurational integral describing intermolecular inter- 
action. The major contributions to the intermolecular 
interactions in water are hydrogen bonding, dipole-dipole 
and quadrupole interactions. The configurational inte- 
gral Z c accounts these interactions. The intermolecu- 
lar interactions in water are stronger than those between 
a water molecule and a field up to the field's strength 
of 2 • 10 9 Vm~ 1A It is assumed that spatial arrange- 
ments of the molecules and density of intermolecular hy- 
drogen bonds are weakly affected by electric field be- 
low 10 9 VmT 1 . Thus the variations of the configura- 
tional integral Z Cl and the translational contribution q tr 
to the partition function with the external field strength 
are neglected. Translational, internal and rotational de- 
grees of freedom contribute to the partition function Q 
of the molecule. The internal degrees of freedom (vi- 
brational and electronic) are only weakly dependent on 
external electric field. In contrast, the molecular dipole 
parallel to the C 2 y symmetry axis of water molecule 
will affect the rotational degrees of freedom in an ap- 
plied field. The above-stated assumptions allow separa- 
tion of the field-independent (translational and internal) 
and field-dependent (rotational) contributions to the par- 
tition function: 



Q = qtr qtn m 



(4) 



where qt r and Qi n are the translational and internal 
molecule's partition functions. The rotational function 
qE for the molecular dipole moment fi in the external 
electric field E can be presented as 

sinh(y) fiE 

<ie = , y = -j—^ ( 5 ) 

y kB-L 

The free energy density change with respect to the state 
E = at arbitrary state of interest is found to be 



RgT . 
-FT— Inte 

Vm 



(6) 



where qt r and qi n are assumed to be independent of the 
electric field. The contribution F 2 to the total free en- 
ergy can be found by integrating the density f 2 over the 
volume: F 2 — J f 2 dV . 
v 

Note that the energy densities fy and f 2 are the func- 
tions of the coordinates (r, z) inside the polymer and 
water meniscus because of the spatial non-uniformity of 
electric field created by the tip-substrate geometry. 

The free energy term associated with the contraction 
of the spring that represents the cantilever bending has 
the form 



U 



k s t 



(7) 



As a first approximation, anharmonism of the spring con- 
stant k s is ignored. 

Finally, the free energy of the tip-substrate electro- 
static attraction depends on the details of the electric 
potential (and the field) inside and outside the film. This 
was determined^ through the exact analytical solution of 
the corresponding Laplace's equation using the method 
of images. 

The energy of the tip-sample electrostatic attraction, 
F ts , was calculated using the image charge construction 
from Ref. [f| . The result has the following form: 



ire V 2 R 2 



E 



^ (n + l)a + h 

n—0 x / / 



(8) 

Here, V is the bias voltage, a is the polymer film thick- 
ness, h = t + R, t is the tip-surface distance, R is the tip 
radius, and rj is given by 



n 



C-2 



ei 



ei + £2 



(9) 



The dielectric constant of water, e 2 , is, strictly speaking, 
a function of the electric field because dielectric satu- 
ration effects become significant at the field strength of 
our interest. For small field's magnitude the water's di- 
electric response is linear with e = 80. The dielectric 
saturation effect becomes substantial as the field grows 
stronger than 10 7 — 10 s VmT 1 resulting in decrease of 
the dielectric constant with the field down to value of 2. 
The molecular dynamics simulations, in particular, pro- 
vide insight into dielectric saturation phenomena. 6 We 
treat t 2 as a constant in the calculations of electric field 
distribution and free energy F 2 . In order to investigate 
the real system behavior, an effective value of e 2 is varied 
between the low and high field limits. 

We use the following dimensionless spatial variables: 
p = r/R, ( = z/R, and the dimensionless electric field 
strength ?/>(/?, £) = ER/V. The following dimensionless 
parameters are defined: 



k s R 



V 2 ' 



w 



RgTR 2 
eoeiV 2 V m ' 
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Rk B T 



(10) 
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After introducing Uq — cq€iV 2 R to scale energy, we fi- 
nally arrive with the formula for the total dimensionless 
free energy function expressed in the units Uq- 



F(t) F 1 F 2 1 o 
— — = — + — + -Kt 2 
U Q U Uq 2 



^ (11) 
U 



where, 



F\ 



a/R R c /R 


is the free energy of the dielectric, and 

(a+t)/R R c /R 

2,t/» / dQ J p In 

o 



F2 



sinh (mip(p, Q) 



mip(p, C) 



(12) 



dp 

is the free energy of water in the meniscus. Eqs. (jllfl . 
© and l|12|) - (|13[) allow calculation of the free energy as 
a function of bias voltage, dielectric properties, and tip- 
sample separation distance. 
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FIG. 2: Dependence of the dimensionless free energy function 
F/Uo on t/R for (1: solid line) bias voltage — 25 V, ei = 3, 
£2 = 80; (2: dot-dashed line) bias voltage — 25 V, e\ = 3, 
£2 = 50; (3: dotted line) bias voltage — 25 V, ei = 3, £2 = 
20; (4: dashed line) bias voltage - 50 V, ej = 3, e 2 = 20. 
The following parameters used: the dipole moment of water 
molecule in condensed phase p. = 1.98 D; the polymer film 
thickness a — R — 35 nm; the cantilever spring constant 
k s = 0.35 Nm" 1 ; T = 25 °C; R c = R. 



III. RESULTS AND DISCUSSION 



The variations of the function F(t)/Uo with respect to 
tip-surface separation t normalized on tip radius R are 
presented in Figure [5] The following typical (reference) 
values have been selected: V = 25 V, t\ = 3, €2 = 50, 
a = 35 nm, and k s — 0.35 Nm^ 1 to analyze the perfor- 
mance of to in poly(methyl methacrylate) (PMMA) and 
polystyrene (PS) films. Starting from the initial refer- 
ence point at close proximity of the tip and polymer, the 
function decreases as the tip retracts from the surface, 
passes through a minimum when tip-sample separation 
is in the order of 30 nm, and then increases again. The 
equilibrium tip-surface separation, to, corresponds to the 
minimum of the function F(t)/Uo. Non-zero value of io 
implies spontaneous tip retraction from the surface, as- 
sociated with the minimization of the free energy. 

The value of the equilibrium tip-surface separation is 
found to be to/R = 0.92. Selecting R = 35 nm, one 
could estimate U = 5.8 • 10" 16 J (3650 eV). The major 
contributions to the system's free energy are due to the 
free energy of water meniscus F2 (— 0.5-10~ 15 J) and 
the spring contraction energy U (0.18 • 10 -15 J) - each 
playing against the other. The contribution of the elec- 
trostatic term F ts (0.11 • 10 -15 J) is comparable to the 
spring energy, whereas the term associated with polariza- 
tion of polymer dielectric F\ is small (— 0.002 • 10 -15 J), 
and does not play a significant role. Thus, the shift in 
the equilibrium tip-surface separation, to, arises from the 
maximization of the volume of water within the highest 
field region under the tip apex - the initial water menis- 
cus lowers its free energy by expanding vertically. 

The value of to depends on parameters such as bias 
voltage V and dielectric constant of polymer t\ (Fig- 
ure EJ, effective dielectric constant of water E2, polymer 
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FIG. 3: Variation of to with respect to the bias voltage (for 
ei = 3, £2 = 50) and with respect to dielectric constant of 
polymer film ei (for bias voltage — 25 V , and £2 = 50). 
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film thickness a and cantilever spring constant k s . An 
analysis of to has been conducted using dimensionless 
parameters, proportional to the first derivative of to with 
respect to the corresponding quantities, calculated in the 
vicinity of the reference values mentioned above. This 
analysis indicates that to mostly depends on the tip bias 
V and cantilever spring constant k s . The variations of 
to with respect to the dielectric constants, and polymer 
film thickness found to be less pronounced. Therefore, 
the nanostructures are expected to be tunable simply ad- 
justing the bias voltage, the easiest parameter to control 
experimentally. This trend seems to be general for dif- 
ferent parameter's values used in this model. 



distance is comparable to the tip's radius. The mechani- 
cal work to lift the tip is produced by the volume of water 
penetrating the tip-surface junction. Another reason is 
associated with the tip repulsion from the surface in the 
double-layered (water and polymer) system. The sign 
of the force acting on the charge Q immersed inside the 
dielectric with permittivity £2 at a certain distance from 
the dielectric with another permittivity e% depends on 
the sign of 62 — ex: repulsive for the positive sign, and at- 
tractive for the negative one. Thus the peculiar behavior 
of the AFM tip can be attributed to energy transfer from 
electrostatic field increasing the potential energy of AFM 
cantilever. This effect presents an additional opportunity 
for high aspect ratio nanostructures formation. 



IV. CONCLUSIONS 



The analysis of the free energy of the system compris- 
ing an electrically biased AFM tip, water meniscus, and 
polymer film indicates that the equilibrium tip-surface 
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